Thermal expansion within a chain of magnetic colloidal particles 
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We study the thermal expansion ol chains formed by self-assembly of magnetic colloidal particles 
in a magnetic field. Using video-microscopy, complete positional data of all the particles of the 
chains is obtained. By changing the ionic strength of the solution and the applied magnetic field, 
the interaction potential can be tuned. We analyze the thermal expansion of the chain using a 
simple model of a one dimensional anharmonic crystal of finite size. 
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It is well-known that the expansion of a solid with tem- 
perature is associated with the anharmonicity of the in- 
teraction between the particles of the solid. This was 
shown with standard solids like lead, using X-ray absorp- 
tion measurements [l[, and in this reference, the thermal 
expansion is analyzed using a one-dimensional anhar- 
monic oscillator model. This suggests that it would be 
valuable to carry out such a study with a one-dimensional 
crystal. Although no true long range order can exist 
in one-dimensional systems in the thermodynamic limit, 
many systems are one-dimensional phases of finite length. 
A classical example is the mercury salts 0] , but more re- 
cently, carbon nanotubes, phases of gases adsorbed inside 
carbon nanotubes Q, or the B phase of DNA have at- 
tracted considerable interest due to their one dimensional 
character. 

In these one-dimensional crystals, only few experi- 
ments have investigated the thermal expansion. One rea- 
son for that is that these systems do not allow for an easy 
control of the interaction between th e particles (or the 
temperature). With colloidal systems, and in particular 
magnetic colloids, it is easy to obtain low dimensional 
condensed phases in which the interaction between the 
particles can be tuned. This strategy was already used 
in where magnetic colloids were used in a beautiful 
experimental investigation of the 2D melting transition. 
To our knowledge, no comparable experimental study of 
ID colloidal crystals has been carried out. 

In the present letter, video microscopy is used to study 
the thermal expansion in a chain of magnetic colloids. 
At high magnetic field, such a chain can be viewed as a 
quasi-one-dimensional condensed phase of colloidal par- 
ticles. In this system, the interaction potential can be 
determined experimentally from the dependence of the 
average interparticle distance with the applied magnetic 
field [3, . We analyze quantitatively these experiments 
using Monte-Carlo simulations. As in standard solids, we 
find that the expansion as function of temperature is con- 
trolled by the anharmonicity of the interaction potential. 
In magnetic colloids, the temperature is measured by the 
inverse of a dimensionless number, A which measures the 
strength of the dipolar interaction with respect to the 



thermal energy. The anharmonicity of the interaction 
potential is varied experimentally by changing the ionic 
strength of the solution, which controls the electrostatic 
repulsive part of the potential. 

In this paper, we describe a colloidal chain as a ID lat- 
tice of N colloidal particles of diameter d with a nearest- 
neighbor interaction potential u(r), where r is the spac- 
ing between two neighboring particles in the chain. Let 
us recall some well-known general features of this system 
before coming to the specific case of magnetic colloids. 
Calculations are conveniently carried out in the ensemble 
where the temperature T and the volume L (which cor- 
responds to the length of the chain in ID) are constant. 
The pressure P (which corresponds to a force in ID) is 
the Lagrange multiplier associated with the constraint of 
constant volume. Thanks to the nearest-neighbors as- 
sumption, the partition function of this problem Z can 
be calculated exactly 0]. It is Z = ^^+7A^, where 
/?= l/fcsT, 



(1) 



and A is the de Broglie length. Thus, the enthalpy G 
equals A^fcBTlog(A/^) (neglecting 1/7V corrections). At 
low temperatures, this system behaves as a solid, which 
suggests to expand the potential u{r) about its minimum 
obtained at rmm- When the function v[r,rmin) defined 
as u{r) — u{rmin), is symmetric with respect to rmin (by 
that we mean precisely w(r, r^m) = w(2r,„„ - r,rmin)), 
we find that there is no thermal expansion (r) ~ rmin- 
At finite temperature, no such symmetry is present: the 
average of the position (r) differs from the most prob- 
able value of the position rmin, which leads in general 
to thermal expansion. This can be seen using a Taylor 
expansion of the form 



v{r,rmin) = ]^kx^ + gx^ , 



(2) 



limited to cubic order with x = r — rmin ■ After calculat- 
ing ^ by the saddle point approximation, one obtains the 
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length L of the chain as a derivative of the enthalpy 
. V P 3gP^ SgkeT 



L 



dP 
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(3) 



The fact that L scales with TV indicates that the expan- 
sion is uniform within the chain. Since we are interested 
in free chains, the pressure P = 0, and the average ex- 
pansion per particle is [J, Q 



(3:) = ^ - rr, 



SgkBT 



(4) 



This equation confirms that the expansion of a solid as a 
function of temperature requires anharmonic terms such 
as the cubic term proportional to 17. It is valid at the 
lowest (linear) order in g or in T. At the same order of 
perturbation, {{x — (x))'^) ~ ksT/k, which is the same 
expression as that obtained for an harmonic potential 
using the equipartition theorem. 

In the specific case of our experimental system of mag- 
netic colloidal chains, the dipolar magnetic part of the 
interaction potential is attractive and of the form 



(5) 



where A ~ jSm? iiij/2i:(P , m is the dipolc strength and /io 
is the permeability of vacuum 0, [l0|- In experiments, A 
is in the range of 10-1000. The non-magnetic part of the 
pair potential is repulsive and can be described as the 
sum of an hard sphere potential Uijs(r) (which is zero 
for r > d and infinite for r = d), and a short-ranged 
electrostatic potential Uei{r). In this paper, we are only 
interested in the regime where the Debyc length is 
much smaller than the particle size. In this case 



Uei{r) =7re*^dln 1 



e 



(6) 



in terms of the zeta potential 5*0, which is assumed to 
be uniform and sufficiently low for the Dcbye theory to 
hold The total pair potential u{r) is 



u{r) = UHs{r) + Uei{r) + Um,{r). 



(7) 



The function u{r) is shown as a function of r on figure 
[1] for different values of A. An important feature of this 
potential is that it is asymmetric for low values of A. 

We have carried out Monte-Carlo simulations of sin- 
gle chains, which give direct access to the thermal av- 
erage (r). Each chain contains 10 to 20 particles in- 
teracting with the potential u{r). We have compared 
results of such simulations to the low temperature ex- 
pansion of Eq. [31 As expected, the low temperature ap- 
proximation correctly describes the simulations with the 
exact potential in the limit where the effective tempera- 
ture 1/A is low. This is shown in Fig. [2] for a particular 
choice of the ionic strength and zeta potential. In the 
range of magnetic field considered and for these condi- 
tions, the low temperature expansion leads to a linear 
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FIG. 1: Total interaction potential it(r) in units of ksT as 
function of the distance r between the centers of two neigh- 
boring colloidal particles. The curves correspond to different 
values of A = 10, 20 or 50 as shown. The black dot represents 
the average position (r) , which is different from the minimum 
of the potential shown as a vertical stroke. The particle diam- 
eter is d = 1150nm, the zeta potential is 'i'o = — 35mV, and 
the ionic strength is 5mM, corresponding to a Debye length 
of about 4nm. 



dependance of (x) as function of the effective tempera- 
ture with (x) ~ 224. 9/A [nm]. Note that error bars on 
the simulation data of figures [5] and [3] were determined 
by the method of Ref. [12l |. they arc very small every- 
where except for the smallest values of A near the limit 
of stability of the chains. 

From a practical point of view, an attractive feature of 
this experiment is that it is possible to use it to measure 
a force acting between neighboring particles. More pre- 
cisely, this can be done by defining the force /((r)), from 
the magnetic part of the potential as 



dum{r) 



dr 



r={r) 



iXkeTd^ 



(8) 



where (r) is evaluated at the middle of the chain. 

This choice of the middle of the chain is made to min- 
imize finite size effects, which we found from our nu- 
merical study have only a small effect on the evalua- 
tion of forces. We have also found that forces calcu- 
lated assuming nearest neighbor interactions differ from 
the case where all neighbors are included [l^, only at 
high magnetic field. In other words, taking into account 
long range dipolar interactions does not change signifi- 
cantly the force versus distance curves and thus is not 
essential for analyzing the thermal expansion effect, as 
we have checked numerically. Furthermore, by compar- 
ing the energy of the finite chain to that of the infinite 
chain, we have found that the force versus distance curves 
depend only weakly on the chain length, thus confirming 
the Zhang and Widom model [3], on which Rcf. is 
based. 
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FIG. 2: Average expansion {x) versus effective temperature 
1/A, the empty circles are the simulation data points with the 
full potential and the filled circles corresponds to the approxi- 
mation of Eq.Uwith a cubic potential. The particle diameter, 
zeta potential and ionic strength are the same as in Fig. [1] 
The solid line is the curve {x) = 224. 9/A which fits well the 
points described by Eq. |3]in this range of magnetic field. 



The force f{{r)) = f{r) defined in this way is equal 
to the derivative of the electrostatic part of the potential 
fei{f) = —duei{r)/dr only when 1/A = 0, i.e. when 
(x) = 0. Thus a deviation of /(r) away from fei{r) can 
be attributed to thermal fluctuations. 

As can be seen in fig.[ll at low values of A, the potential 
u(r) is very steep for r < where it can be approx- 
imated by its short-ranged part Uei{r) + UHs(r), while 
u(r) ~ Umir) for r > rmin- The most asymmetric situa- 
tion corresponds therefore to the limit where the Debye 
length goes to zero. In this case, the magnetic part of 
the potential is only balanced by the hard sphere repul- 
sion of the colloids. This case is more simple than the 
previous case because now the position of the minimum 
of the potential rmin is always d irrespective of A. An 
expansion of the potential about this minimum a,t r — d 
leads to 



1, 

v(x) ^ px -\ — kx 



(9) 



with p — u'm{d) — i\kBT/d and k ~ u'^{d) = 
— 12XkBT/d'^. Now by treating the term kx^/2 as a per- 
turbation with respect to the first term px, one obtains 
to lowest order in T: (x) ~ ksT/p, and thus 



(r) ~ d 1 + 



1 

3A 



(10) 



This approximation holds for A » 2/3, so essentially for 
all values of the magnetic field considered here. At the 
same order in perturbation, {{x — {x))"^) — {kBT/p)^, 
which means that in this case the RMS fluctuations of x 
are of the same order than their average. For hard-sphere 
potentials, Eq. [TOlcan be inverted to obtain A((r)) which 
leads to the following simple expression for the force 



(11) 



In the limit T = 0, we recover that /((r)) 
3fiom^ /2iTd'^. The force /((r)) shown in Fig. [31 has a 
vertical part near (r) ~ d and a slower power-law decay 
at larger distances. The vertical part of the curve typical 
of hard sphere behavior has been observed in Ref. [lB| 
but the slower decay at larger distances has not been an- 
alyzed there although it is expected to be present in these 
experiments. 

Let us discuss measurements of the force versus dis- 
tance [5|. We use Dynal-Invitrogen magnetic beads of 
diameter d ~ 1150nm, as observed by Dynamic Light 
Scattering. The polydispersity in size is estimated to be 
in the range of 1%. We vary the interaction potential by 
changing the ionic strength, and we measure the mean 
particle distance using digital video- microscopy [l^ . The 
mean distance is calculated by averaging the distance be- 
tween particles within the chain and over the time of the 
experiment. The force versus the mean distance between 
bead surfaces is shown in Fig. [S] Each set of experi- 
ment with a given ionic strength is performed with the 
same chain in order to reduce errors due polydispersity 
in particle size. To determine the electrostatic poten- 
tial, we have fitted the quasi-linear part of the curve for 
the zeta potential "ifo, and we have found ^'o — — 35mV. 
When the Debye length is large, the expansion (x) is al- 
ways negligible in the range of force and distance which 
can be resolved. As a consequence, we observe only the 
low temperature regime with a quasi-linear force (in log 
scale) versus distance profile, characteristic of the electro- 
static part of the interaction. For smaller Debye length, 
thermal fluctuations make the expansion (x) observable. 
A very good agreement is found between the simulations 
and the experiments except for the smallest Debye length 
in the region where the deviation from the linear profile 
occurs. We believe that this discrepancy is not due to a 
failure of the electrostatic description, which could only 
happen at much higher values of the surface potential ^'o 
than observed here. Instead we attribute the discrepancy 
to an additional repulsive force, which is due to a polymer 
layer grafted on the particles 17[. This force is indeed 
important only when the particles are sufficiently close 
but irrelevant otherwise. When this additional contribu- 
tion is taken into account, a good agreement between the 
experiments and the theory can be obtained for all the 
ionic strengths studied here. 

In Ref. 0, a deviation in the force versus distance 
measurements (with respect to a linear profile) was re- 
ported at high ionic strength and for low field, which 
we attribute to the thermal expansion effect discussed in 
this paper. The particles were a factor 5 smaller (there 
d ~ 190nm) than in the present study. In fact, for smaller 
particles, a higher ionic strength is necessary to observe 
the deviation of Fig. [31 The importance of the parti- 
cle size or of the Debye length can be determined by 
evaluating the position of the point where a deviation 
starts to become observable. If we assume that this hap- 
pens when (x) ~ 0.1 x (r), using our low temperature 
approximation, we find that this point corresponds to 
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FIG. 3; Force versus the mean distance between the bead sur- 
faces (r) — d for different ionic strengths. Symbols correspond 
to experimental data points, from left to right, ▲ 5mM, • 
ImM, ^ 0.25mM and ■ 0.07mM. The straight lines represent 
the electrostatic part of the force, the dot-dashed lines repre- 
sent Monte Carlo simulations, and the dashed line on the left 
represent the expected force for a hard sphere potential. 

[F = 0.7pN,/i = 30nm) for the curve at c = 5mM and 
{F = 0.17pN,/i = 131nm) for the curve at c = 0.25niM 
in agreement with Fig. [31 For lower ionic strength and 
lower forces, the effect is too weak to be observable. 

As the strength of the dipolar interaction is reduced, 
we observe both in the experiments and in the simula- 
tions that thermal fluctuations cause chains to fragment. 
Independent of fragmentation, disordering begins at the 
ends [3, where the inter-particle distance is larger due 
to finite size effects. This effect is similar to edge and 
surface melting in 2D and 3D solids. We have stud- 
ied numerically the ratio of the standard deviation of 
the inter-particle distance divided by the average inter- 



particle distance, L^. = \J {x^) / {r), which is similar to 
the Lindemann parameter used to detect the approach 
to melting in 3D systems. At large A, the potential is 
harmonic and is small, whereas for small A, the po- 
tential becomes anharmonic and Lg increases. We find 
that Lf, raises very steeply as it approaches a value of 
the order of 0.1, which corresponds approximatively to 
melting in 3D systems. The threshold value on A where 
this happens depends very much on dimension, it is much 
lower in 2D as in ID for instance. 

In order to investigate the 2D case, we have developed 
a 2D extension of our Monte Carlo simulation, in which 
the location of the magnetic dipoles of the chain form 
a 2D polymer. In addition, we have assumed that the 
orientations of these dipoles remain frozen in the direc- 
tion of the applied magnetic field. We have found that 
the force versus distance curves shown in Fig. [3] do not 
differ in the 2D case as compared with the ID case, in 
the range of magnetic field considered here. A difference 
between the predictions of the 2D and ID model only 
arise at very small forces and large distances. In the 2D 
case, transverse modes of fluctuations of the chain arise, 
which are soft bending modes [T^]. Such soft modes do 
not significantly alter the mean distance, which is repre- 
sented in Fig. [21 This explains the robustness of the ID 
model. The region where thermal fluctuations are very 
large requires further modelling, since there the 2D or 3D 
character of the problem is important. 

To summarize, we have analyzed the thermal expan- 
sion of a chain of magnetic particles. We have shown 
that this effect is responsible for a deviation with respect 
to the quasi-linear force profile observed at low effective 
temperatures. This effect can be described using a simple 
ID model, and is relevant to the determination of force 
versus distance using magnetic colloids. 
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